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ABSTRACT 



Context. Exploring the structure and dynamics of cold starless clouds is necessary to understand the different steps leading to the 
formation of protostars. Because clouds evolve slowly, many of them must be studied in detail to pick up different moments of a 
cloud's lifetime. 

Aims. We study here a fragment of the long filament L1506 in the Taurus region which we name L1506C, a core with interesting 
dust properties which have been evidenced with the PRONAOS balloon-borne telescope. 

Methods. To trace the mass content of L1506C and its kinematics, we mapped the dust emission, and the line emission of two key 
species, C'**0 and N2H+. ''CO and C"0 were also observed. We model the species emission using ID Monte Carlo models. 
Results. This cloud is reminiscent of LI 498 but also shows peculiar features: i) a large envelope traced solely by ''CO holding a 
much smaller core with a strong C'*0 depletion in its center despite a low maximum opacity (Av~ 20 mag), ii) extremely narrow 
C'^O lines indicating a low, non-measurable turbulence, iii) contraction traced by C'**0 itself (plus rotation), iv) unexpectedly, the 
kinematical signature from the external envelope is opposite to the core one: the "CO and C"*0 velocity gradients have opposite 
directions and the C'^O line profile is blue peaked on the contrary to the "CO one which is red peaked. The core is large (r = 3 x 10* 
A.U.) and not very dense (n(H2) < 5 x 10* cm"', possibly less). This core is therefore not prestellar yet. 

Conclusions. All these facts suggest that the core is kinematically detached from its envelope and in the process of forming a prestellar 
core. This is the first time that the dynamical formation of a prestellar core is witnessed. The extremely low turbulence could be the 
reason for the strong depletion of this core despite its relatively low density and opacity in contrast with undepleted cores such as 
L1521E which shows a turbulence at least 4 times as high. 

Key words. Stars: formation - ISM: clouds - ISM: abundances - ISM: molecules - ISM: Structure - ISM: individual: L1506 



1. Introduction 

On their way to form stars, clouds go through a contraction phase 
' leading to the format ion of prestellar cores (n > 1 x 10^ cm"^, 
: iKeto&Caseliill2008h which are usually identified by their lack 
' of internal heating sources, their large CO depletion, and their 
, strong NH3 and N2H^ emission lines for which turbulence is of- 
ten subthermal. These cores subsequently collapse to form pro- 
tostars. Many details pertaining to these two condensation steps 
are still obscure and we don't understand what conditions are 
necessary to produce stars of different masses, high or low, what 
conditions are necessary to start the contraction of the cloud or 
the collapse of the prestellar core, nor the role of the magnetic 
field in supporting the clouds against collapse, etc... Obviously 
as the evolution time is long, one has to observe many clouds as a 
substitute to tracking any single cloud in the process of forming 
a star For each cloud, physical and chemical properties should 
be studied in detail to assess all important parameters such as 
temperature, density, kinematics, and their respective gradients, 
if any, and also chemical abundances by volume. By accumu- 
lating this type of studies, we can hope to disentangle depletion 
phenomenon versus cloud age or density, contraction versus os- 



* based on observations made with the IRAM-30m. IRAM is sup- 
ported by INSU/CNRS (France), MPG (Germany), and ION (Spain). 



dilation, etc. and describe how prestellar cores and protostars 
form. 

Many prestellar cores have been identified and some have 
already been studied in detail, including B68, LI 83, L1498, 
L1544, etc. Among these objects a few cases have been identi- 
fied as early prestellar cores if their density is already greater 
than 1 X 10^ cm"^, or simple cores probably about to form 
prestellar cor es if not. T he first type includes L1495B, L1521B, 
and L1521E ( Hirota et al. 2002 . 2004 ; Tafalla & Santiag o 2004|) 
which, despite a high density, show no sign of depletion (strong 
CO lines, weak or absent NH3 and N2H^ lines) and are con- 
sidered as being both physically and chemically young while 
the second type is represented by two cases presently, L1498 
jTafalla et al.ll2004l;IShirlev et alj|2005l) and the verv d ense glob- 
ule ' G2' in the Coalsack comple x TLada et al.ll2004l) . L1498 in 
the study bv ' Shirley et alj (12005') exhibits a low density (1-3 x 
10'* cm"^) and molecula r depletion ( al ready reported by several 
other authors, especiallv lWillacv et ar][T998l who, by studying 
L1498, reported the first C'^O depletion in a prestellar core by 
comparison with a dust map). The G2 gas content has not been 
studied in detail yet and only its peculiar, ring-like dust shape 
indicates its unstable nature despite its low density. Thi s core is 
expected to turn into a prestellar core (iLada et al.ll2004 . In this 



second category, which is considered as physically young but 



2 



L. Pagani et al.: L1506 : a prestellar core in the making 



chemically evolved, we present a newcomer, a fraction of the 
LI 506 filament. 

The cloud we present in this paper, LI 506, ha s attracted at- 
tention during a large-scale '^CO survey (Aberge retanil994 
because of its low (though not unique) DIRBlJ] Ii40/jm/l240//m 
and IRAS0 leo/jm/Iioo/jm colour ratios to which the '^CO maps 
were compared. Its dust properties were studied later on using 
the PRONAO^ balloon-borne experiment contin uum emission 
data i n the submillimeter range (200 - 600 fim) jStepnik et al.l 
I2OO3 I). albeit with a limited angular resolution (2-3.5')- A sig- 
nificant change of the dust properties has been revealed going 
from the diff'use to the dense part of the cloud : no emission 
from transiently heated small particles, and a strong enhance- 
ment (by a factor of 3.4) of the submillimeter emissivity in the 
latter. This was interpreted as the signature of dust coagulation 
leading to the formation of fluffy aggregates (made of a mix- 
ture of very small and big grains). In order to understand this 
evolution and investigate what are the physical conditions asso- 
ciated with efficient aggregate formation (ice mantles, level of 
turbulence, densities,), we are now concentrating on the study of 
the gas properties. It is presently clear that depletion plays a key 
role in the gas phase of dark clouds and basically two species are 
needed to map the clouds, namely C"*0 and NaH^. NH3 is also 
a good tracer but N2H^ has some advantages over NH3 among 
which the possibility to use the same telescope as for C"^0, 
giving a similar resolution and beam corr ection and also supe - 
rior diagnostic capabilities as discussed bv lPagani et al.l (l2007h . 
By observing these species, both the undepleted (outer) and the 
depleted (inner) cores can be traced and their kinematics and 
physical properties can be modeled. With '^^CO , we will also be 
able to trace a fraction of the extended envelope. We present and 
model such observations in this paper and discuss the properties 
of this cloud. 



2. Observations 

All observations were performed using the IRAM-30m tele- 
scopjB- The dark nebula L1506 is an elongated object 
(IOnisliietal.llT996t iNercessian etaDflgSSh and the SIMBAE0 
reference position (oraooo = 4''18'"31.r ^2000 = -h25°19'25") is 
outside the part of the cloud we present here which is cen- 
tered on the P RONAOS emission p eak (a2000 = 4''18'"50' (520oo 
= +25°19'15". IStepniket aT]|2003l) wh ich we propose to n ame 
L1506C (A and B are already defined, 'Lee & Mvers"l9 9l. Its 
distan ce is estimated to be 140 pc (Elias 1978; Kenvo net alj 
11994 . 

In September 2004, we mapped the filament dust emission at 
1 .2 mm with the MAMBC0 II bolometer. Sky dips, pointings and 
calibration source observations were regularly performed, and 
data reduction was executed using IRA M proprietary software. 
The cold filament studied bv lStepniket al.l (1200 3) is clearly vis- 
ible at the centre (Fig.[T]). Spectroscopic observations were per- 
formed over several runs: 29 & 30 November 2003, 22 & 23 May 
2004, 28 July 2004 and 18 July 2008. In all cases, the VESPAQ 
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Fig. 1. Dust emission at 1.2 mm as measured by MAMBO II. 
Resolution has been degraded to 30" for a better signal-to-noise 
ratio. Levels are spaced by 0.5 MJy sr"' from 0.5 to 4 MJy sr"' . 
The dashed line represents the cut observed with PRONAOS, 
and the white section the fraction of the cut observed with the 
IRAM-30m. The original cut direction was chosen to be per- 
pendicular to the FIR filament (colour image on-line). The (0,0) 
position corresponds to q;2000 = 4''18'"50' ^2000 = -i-25°19'15". 



autocorrelator was used with frequency sampling varying from 
6 to 20 kHz for 3 mm lines and 40 kHz for 1.3 mm lines (Table 
[TJ. All observations were done in frequency switching mode. 
Receivers A&B were used except in July 2008 where the 1 mm 
9 pixel cameras HERA0 1&2 were used in the On-the-Fly mode, 
still using frequency switching. The C"^0 and '^CO (J:l-0) lines 
were observed simultaneously by tuning the receiver halfway be- 
tween the two frequencies and off'setting the VESPA autocorrela- 
tor subwindows by +210 MHz in front of the two lines. Pointing 
was regularly checked and found to be stable within 3". 



Table 1. Telescope settings 



Line 



Frequency 
(MHz) 



Beam size 
(") 



Sampling 
(kHz) m s" 



N2D+ (J 


1-0) 


77 109.616" 


32 


6 


23 


N2H+ (J 


1-0) 


93 173.764" 


27 


10 


32 


c"*o (J 


1-0) 


109 782.176 


23 


10 


27 


"CO (J 


1-0) 


110 201.360 


23 


20 


54 


C"0 (J 


1-0) 


112 358.982" 


23 


10 


27 


C"*0(J 


2-1) 


219 560.358 


12 


40 


54 



" New frequency for the m ain 1^^/ : I23 
from lPagani et alj | |2009 
Jf:17/2-05/2 transition from lKlapper et alj ( |2003|) 



O12 hyperfine component 



HEterodyne Receiver Array, see IRAM website 
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Fig. 2. Left: low res olution '^CO (J: 1-0) integrated intensity map of the L1506 Taurus cloud filament (original data from 
iNercessian et al. l [T988h . Contour levels are 0.5 to 4.5 Kkms ' by steps of 1 Kkms ^ The POM-1 2.5-m dish has a 4.5' reso- 
lution as shown in the bottom right corner The square indicates the surface mapped in C'^O (J:2-l). The offsets are with respect to 
L1506C coordinates (see text). Right: the smoothed C''^0 (J:2-l) integrated intensity map. Superposed to the map is the strip path 
observed with PRONAOS (Step nik et al.l l2()03) along which we observed '^CO, C'**0, and NaH^ with a good signal to noise ratio 
(see Fig. O. The crosses indicate the observed points and the 3 white boxes (superposed to 3 crosses) the points observed also in 
C'^O (see Fig.|5]l. The degraded 30-m beam size is represented. Levels are from 0.1 to 0.7 Kkms"' by steps of 0.1 Kkms '. The 
(0,0) position corresponds to Qf2oo() = 4''18'"50' ^2000 = +25°19'15". 
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Fig. 3. Spectra taken along the 
PRONAOS cut (see Fig. The 
upper row shows the '^CO (J:l- 
0) spectra. The middle row shows 
the C'^O (J: 1-0) spectra. The lower 
row shows the N2H^ (J: 1-0) spec- 
tra. The velocity scale is 6.5 to 8 
kms ' for '^CO and C'^O, and 5.5 
to 9 kms ' for N2H+. Offsets in- 
crease from North- West to South- 
East. 



3. Results and analysis 

Figure [2] shows a large-scale '^CO (J: 1-0) map obtained with 
POM-lQ (a 2.5 m dish telescope with 4.5' resolution and 100 

' Petite Operation Millimetrique 



kHz filters - 0.27 kms ' - which was in service at Bordeau x 
Observatory. Original data taken from INercessian et al.lll988l) . 
and the region of interest is delineated by a square. This map 
traces the large-scale envelope of the cloud. On the right hand- 
side, the C"^0 (J:2-l) integrated intensity map obtained on-the- 
fly with the 30-m over L1506C is shown. We degraded the res- 



4 



L. Pagani et al.: L1506 : a prestellar core in the making 



C O 
N,H* 




offset (orcsec) 

Fig. 4. Integrated intensity of C'**0 and N2H+ (J: 1-0) com- 
pared to the dust emission as measured by MAMBO along the 
PRONAOS cut. MAMBO data are smoothed to 40" here. The 
core is clearly not symmetrical and C'^O depletion is obvious. 

olution to 24" to improve the signal-to-noise ratio. A large hole 
in emission is visible in the centre of the map. This is due to 
CO depletion, as is now well established. The PRONAOS strip 
centered at the 100 fj.m emission peak came across the most de- 
pleted part of the filament. This depletion is also spectacular in 
terms of spectrum peak temperature (Fig. [3]) or integrated inten- 
sity (Fig. nil along the PRONAOS strip. In these figures, the drop 
of C"^0 intensity is correlated with the appearance of N2H^, an- 
other sure sign of CO depletion as CO readily destroys N 2H^ 
( lBergin&Langeijll997l: ICasellil I2OO2I: iPagani et alj|2005h . To 
check that the C^^O lines are not dropping for opacity reasons 
in the centre of the cloud, we observed 3 positions in C'^O (J:l- 
0) . We chose two positions close to the C"^0 peaks on each 
side of the strip (+75") and one in the middle, close to the lo- 
cal minimum of the C'**0 emission (Fig.|5]l. These 3 positions 
are indicated on Fig. |2]by 3 small white boxes superposed to 
3 of the crosses along the PRONAOS strip. Several interesting 
features are visible: 

1. the (0, 0) and (+75", 0) offset positions (along the cut) dis- 
play extremely narrow C'^O and C'**0 lines (FWHM = 0.16 
to 0.26 km s ' where FWHM stands for full width half max- 
imum). 

2. the existence of a velocity gradient across the strip which 
possibly traces rotation in the C'^'O data. 

3. the existence of another velocity gradient for '^CO of oppo- 
site direction. 

4. the C'^O (J: 1-0) splitting towards the reference position, 
which is not so much visible in the (J:2-l) line but is clearly 
so in the C'^O (J: 1-0) line, as the dotted lines reveal in Fig. 

s 

As most of the properties of this cloud change gradually with 
radius, it is difficult to define what is the core and what is the 
envelope. One parameter however changes abruptly and this is 
the velocity field. We will therefore define the envelope as the 
region where '^CO is revealing a different velocity field from 
the other species while C''^0 and NaH^ both trace the core, the 
outer, undepleted part for the former and the inner depleted one 
for the latter 

As far as we know, these CO narrow lines are the narrow- 
est C'**0 and C'^O lines reported so far, comparable to the 
very narrow features seen in NH3 and N2H^ in sources like 



L183 dPagani et al.l2007h . The main difference is that these C"*0 
and C'^O lines are emitted outside the depleted core, in a re- 
gion where turbulence is usually still high. The lines are narrow 
enough so that the C'^O Jf:17/2-05/2 and Jf:13/2-05/2 hyperfine 
components which are usually blended, are completely separated 
here. 

Figure |6] confirms the extraordinary narrowness of the lines 
towards the centre of the cloud. These C"^0 spectra have been 
taken at (0,-25") in the frame of the PRONAOS cut, that is (- 
21.6", -12.6") in the Equatorial J2000 frame. The (J:l-0) spec- 
trum shows clearly two components which we have fitted with 
2 independent gaussians. We have also fitted the (J:2-l) spec- 
trum with two unconstrained gaussians. The velocity of each of 
the two components agree for both lines within the uncertainty 
of the fit (a difference of 10 ms ' for an uncertainty of 30 and 
40 ms~' for the left and right components, respectively). As for 
C' ■'O in Fig.|5] the splitting between the two components is -200 
ms"' and the width of the individual components is 150 (right) 
and 180 (left) m s ' . For a standard gas kinetic temperature of 10 
K, for which the optically thin thermal linewidth of CO is ~0.12 
km s ', this sets the turbulent width contribution to the linewidth 
to <90ms-'. 

If the C''^0 and C'^O splitting was due to self-absorption, 
it would be exactly the reverse, the C"^0 (J:2-l) line would 
have had the deepest self-absorption feature and normally no 
such feature would have been seen in C'^O since, due to its 
low abundance and to its hyperfine structure, it is very difficult 
to make it optically thick. These symmetric features cannot be 
due to rotational Doppler shifting because for a narrow beam 
aiming at the rotational axis, most of the material picked up 
by the beam is moving at right angle to the line of sight hav- 
ing therefore no contribution to the velocity component along it. 
Only material very near the rotational axis would contribute to 
the Doppler shifted emission. As the inner core is depleted in 
CO, such a contribution near the axis is not possible and there- 
fore rotation cannot explain this splitting. To explain this double 
peak feature, the only possibility left is that the CO envelope 
is radially moving. The direction of movement (expansion or 
contraction) can be told from the asymmetry of the C'^O (J:l- 
0) and (J:2-l) lines. Here, these C"*0 lines are stronger in the 
blue component than in the red component. As we will confirm 
with the Monte-Carlo model (Sect. 4), this can only be explained 
by contraction and not expansion despite the absence of heavy 
self-abs orption which is us ually invoked to justify blue peaked 
profiles jMvers et al.lll996h . The lines being almost independent 
due to their narrowness, the blue and red line contributions are 
emitted from only one side of the cloud each - front or back de- 
pending on the movement, expansion or contraction - and have 
no or very little interaction. The line which comes from the rear 
part is originating in the backside low density envelope and gets 
contribution from a material getting denser as the photons travel 
towards the observer, until they reach the depleted zone and then 
cross the rest of the cloud without interaction. For a constant ki- 
netic temperature, this means that the excitation temperature of 
the contributing material is increasing (or remaining constant if 
thermalized) towards the observer and the successive layers add 
their contribution to the emergent signal. Conversely, the line 
formed on the near side of the cloud starts from the highest den- 
sity and excitation temperature layer, just outside the depleted 
zone and gets contributions from lesser and lesser dense mate- 
rial with therefore a decreasing excitation temperature. When the 
excitation temperature decreases, net absorption of photons can 
happen if opacity is not null. Indeed, the C '^O lines are not really 
optically thin (opacities are in the range of 0.1 to 0.5) and this 
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Fig. 5. Spectra of '^CO (dotted lines), and C''^0 taken at 3 

different positions along the PRONAOS strip: in the middle and 
close to each side peak C'**0 emission. The full width at half 
maximum is indicated for some remarquably narrow lines. Only 
the C'^O Jf:17/2-05/2 and if-^j/i-^s/i are displayed. The dashed 
line marks the systemic velocity (with an offset of +0.534 km s"' 
for the Jf:13/2-05/2 C'^O transition) while the dotted lines for the 
C'^O transition show the symmetrical velocity displacement of 
each C^^O hyperfine component in the (0,0) direction (see text). 
Lines are artificially offset vertically and the central offset C'^O 
line is multiplied by a factor of 2. The C'^O (J:2-l) lines are 
shown at their original 11" resolution. 



(0,-25") 




6 7 8 9 

LSR Velocity (km s ) 

Fig. 6. C'^O (J: 1-0) (top) and (J:2-l) (bottom) spectra (his- 
togram plot) taken 25" away from the PRONAOS cut. The total 
Gaussian fit (full line) for both transitions is shown as well as 
individual components. The components are separated by 200 
m s"' and their width is 150 (right) and 180 (left) m s"' 



is enough to slightly differentiate between the line which comes 
from the rear part of the cloud which has undergone no attenua- 
tion after the highest excitation temperature layer emission and 
the line which comes from the front part and has suffered from 
a small absorption due to the outer, less excited layers. The ef- 
fect is clearly seen towards the (0,0) position : the C'**0 (J: 1-0) 
line is 20% stronger on the blue side (1.81 and 1 .48 K peak tem- 
perature for the blue and red components, a difference greater 
than 10 0-, with a ^ 21 mK), and the C'^O (J:2-l) Une which 
due to higher opacity and higher sensitivity to density condi- 
tions has lost its red peak (Fig.|5]). The red component is there- 
fore the one emitted from the front part of the cloud, the blue 
component from the rear part and this indicates that the cloud is 
contracting rather than expandin^3 The amount of contraction 
of the C'^O core is indicated by the two dotted lines in Fig.|5] 
namely ^ +100 ms"'. Therefore, the C'^0 outer core both ro- 
tates (from tlie (+75",0) symmetrical velocity offsets seen in 
Fig. m as mentioned above) and contracts. To our knowledge, 
this is the first time that core contraction is observed based on 
C'**0 and C'^O lines. 

The N2H^ hyperfine structure lines do not have such a 
good signal-to-noise ratio and do not show a structured profile. 
However, their displacement in velocity along the cut is similar 
to that of C'**0 as we can see in Fig.]?] 

The '^CO (J: 1-0) line which was observed simultaneously 
with the C'**0 (J: 1-0) line shows a very different behaviour 
(Figs. [3] 12] &|7]i- Its peak intensity is red-shifted and its veloc- 
ity gradient along the cut is opposite to that of the C"*0 line. 
From the large scale '^'CO map obtained with POM-1, we know 
that the L1506 filament has a width of ~40' in the region around 
L1506C (Fig. ID. This is ~6 times larger than the C'**0 extent. 
As we do not know how the '^CO velocity behaves beyond these 
central 6', we cannot claim that this velocity drift is the signature 
of a rotation of the large scale filament but it definitely indicates 
that it is opposite to the C'**0 rotation (Figs.|5]&|7]l. Similarly, 
the red peak in the '^CO velocity profile is probably a sign of 
gas expansion (Figs.[3]&|5]l, while C'^O is tracing a contraction. 
Therefore the extended envelope is dynamically behaving at the 
opposite of the core. 

Figure |4] shows the dissymmetry of the cloud in both dust 
and gas. It is remarkable that the C'^O lines on each side of 
the cut have the same maximum integrated intensity despite a 
variation of 25% in peak intensity. It indicates that the quantity 
of gas is comparable and that the main reason for the intensity 
change is the line width difference. It must also be noticed that 
the N2H^ integrated intensity profile is strongly asymmetric and 
varies abruptly at the top, two features we cannot reproduce with 
our standard one-dimensional (ID) model (see Sect.|4]i. Its peak 
position is close to but not coincident with the dust peak position. 

We observed the position with the second strongest N2H^ 
integrated intensity (offset +25" along the PRONAOS strip) in 



There exists the alternative possibihty that the blue peak is stronger 
than the red peak because of a difference in column densities between 
front and rear parts of the cloud. If it were the case in such a manner 
that it would hide the cloud expansion, the difference in column density 
between the two components would be ~40% in the optically thin limit 
(and even more with growing opacity). Indeed, as the intensity differ- 
ence is about 20% towards the (0,0) position in (1:1-0), it would 
have to be compensated twice to go from a stronger red component 
(due to expansion combined with differential excitation) to a weaker red 
component (as observed) and the effect would be more pronounced for 
the (1:2-1) component. This possibility seems therefore quite unlikely 
and even if the two sides of the cloud have different column densities, 
contraction remains the best explanation. 
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^^CO (1-0) (1-0) (2-1) N2H^(1-0) 

Observations Model Observations Model Observations Model Observations Mode 
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Fig. 7. Position-velocity plots along the PRONAOS cut. For each line, we show the observations on the left and the model on the 
right with the same contour levels. For N2H^, we reproduce only the 3 inner hyperfine components. The dashed line marks the 
systemic velocity of the cloud. For '^CO and C'**0 observations, the dotted contours indicate interpolated data (data is missing for 
offsets +162.5")- Superposed to the '^CO observations, dotted contour levels from the C'^O (J: 1-0) observation plot are repeated. 
Only levels 1, 2, and 3 K have been drawn for clarity. 



N2D^ (J: 1-0) to check the level of deuteration in the inner core 
(the strongest peak - offset +37.5" - was identified afterwards). 
We found only a very weak signal (Fig. [8]i, the main hyperfine 
component being only 70 mK. A hyperfine fit was performed 
using the CLASS fitting routine MINIMIZE with the HFS op- 
tioiO 



4. Models 

4.1. density profile 

The density radial distribution has been derived from the 
MAMBO map emission using a Wiener linear inversion method 
as described by Dupac & Giard (2002). Considering the elon- 
gated distribution (along the RA axis), we have modelled the fil- 
ament with a cylindrical geometry approximation, with its axis 
in the plan e of the sky and pe rpendicular to the PRONAOS 
strip, as in Stepni k et al.l (l2003l) . We have also taken into ac- 
count a temperature inward decrease from 15 K down to 8K. 
This r ange has been dedu ced from the Av profile in the fila- 
ment (Stepnik et al jj2003[) combined with th e prediction from 



Bernard et al.l d 19921) dust temperature 



iZuc coni et al. (200JJ) or 
models in cold cores. We have discretized the cylinder into iso- 
density and iso-temperature rings, and determined the array A 
such as: 

Nmd) = A(i, i) X riHiii) + b 



or 



where (i,j) are respectively the indexes for the rings and the lines 
of sight, Nh2 the column density, nH2(i) the gas density in the 
ring i, /j the observed intensity along the line of sight j, b is the 
noise level, Bi(T,) is the black-body function evaluated for the 
temperature T of the ring i, m// is the proton mass, - 2.33 



http://www.iram.fr/IRAMFR/GILDAS/ 



is the mean molecular weight of interstellar material in molecu- 
lar clouds and /C1200 /im = 8 x lO^-'cm^^g"' is the mass absorp- 
tion coefficient assumed constant in the core and derived from 
the dust aggregates emissivity model by lOssenkopf & HenningI 
( 1994). The local density nH2 in each ring (i) is deduced from : 
tiH2 - W X I^, where we adopt the optimized inversion matrix : 

W = ('A' xA')-^ x'A' 

(cf Dupac & Giard (2002)), with: 

A' = /^i2oo ^lm X niH x A{i, j) x BAiJi) 

The best optimization has been obtained using a ring thickness 
value of 50", and the result leads to an inner density of 2.4 x 10"* 
cm"-' (in average over a 50" radius core). 



4.2. line modeling 

The density profile derived from the dust data has been used as 
a first guess for our ID line radiative transfer models but we 
allowed the profile to be changed because the raw dust map is 
noisy and its final resolution for the inversion method is only 50" 
compared to the C"*0 (J:l-0) 24" and (J:2-l) 12" and because 
the dust properties are not strongly constrained and can vary, 
therefore making the derived density uncertain by a factor of 2 
typically. 

For all three species, '^CO, C'^O, and N2H+, we fit a 
non-local thermodynamic equilibrium radiative transfer Monte- 
Carlo model and try to reproduce the spectra (Figs. |9l [TOl 
nn & [T2I 1 and the velocity-position diagram (Fig. [T]!. The 
ID M onte-Carlo model was originally developed bv BernesI 
(Il979h . It includes microscopic turbulence and a radial ve- 
locity field. It was subsequently modified to include rotation 
(iPagani & Breart de Boisangeii [l996) and a varian t was devel- 
oped to treat the hyperfine structure of N2H'^ jPagani et al.l 
2007*). The ID cloud parameters are shown in Table |2] and the 
species abundances are traced as a function of radius in Fig. 
[13] Here the model is spherical, with a radius equal to the dust 
cylinder radius and with the rotation axis similar to the cylinder 
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Fig. 8. Spectrum of N2D+ taken at offset +25"along the 
PRONAOS strip. The hyperfine fit is superposed to the data. The 
data is Harming smoothed once. The sampHng is 50 ms ' and 
the rms is 1 3 mK. 

axis of the dust model, i.e. supposed to be aligned with the em- 
bedding filament elongation and therefore perpendicular to the 
PRONAOS strip. The model data to be compared to the obser- 
vations are taken at the equator considered to be a good approx- 
imation of a cylinder slice, as already discussed in Pagani et al. 
(l20Q7h for a similar case. For the '^CO emission, as our high 
spatial and velocity resolution strip only extends upon 6' and be- 
cause the '^CO line is optically thick and much more extended, it 
is difficult to build a simple ID model to reproduce and constrain 
the 30-m '^CO observations (the POM-1 data are too coarse and 
too noisy for this purpose). More generally, the cloud is not sym- 
metrical and therefore our ID Monte-Carlo model cannot repro- 
duce the various differences between positive and negative off- 
sets. However, we have tried to fit globally the observations and 
overall to use a similar cloud description for the N2H^ and C'^O 
line models. Though the '^CO data are too difficult to reproduce 
with the ID Monte-Carlo code, we have attempted to fit them as 
well, at least to get some estimate of the physical conditions in 
the outer parts which could explain the observations. From the 
C"^0 and N2H^ data, we find that the cloud rotational axis could 
also be slightly offset along the PRONAOS cut (about -1-12") 
with respect to our reference position (which was set somewhat 
arbitrarily). The peak N2H^ integrated intensity is even further 
away from the reference position (+37.5") but setting the rota- 
tional axis that far does not fit with the overall velocity profile. 
Despite the fact that the cloud extension on the positive offset 
side is smaller than on the negative offset side, the fit is relatively 
good, especially towards the centre of the cloud. As mentioned 
above, we have tried to somewhat adjust the density pr ofile and 
found a solution which is close to iTafalla et al.l (l2002h formula 
for the N2H+ - C"*0 region : 

«(H2) = ---7 cm (1) 

with no(H2) = 5 X 10"^ cm-\ ro = 1 x lO^A.U. and a = 2.5. 
This profile provides a peak dust emission comparable to the 
MAMBO observation, and the average density of ~4x 10"* cm""* 
in a 50" radius is within a factor of 2 of the peak density derived 
from the MAMBO observations. Apart from the dust properties 
which are not strongly constrained, an equivalently possible ex- 



planation for this factor of < 2 differenc e is that the He + N 2H"^ 
collisional coefficients we presently use dDaniel et al.ll2005h are 
overestimating the true H2 density needed t o colhsionally excit e 
the N2H^ to its J = 1 level as discussed in iPagani et alJ (l2007h . 
From dust and N2H^ measurements, we find that the density in 
this core is rather low, a factor o f 2 to 4 below the threshold pro- 
posed by iKeto & Casellil (l2008h for prestellar cores. Therefore, 
it is not such a core. 

The C0+N2H^ model which we present here is by no means 
unique but represents a family of possible solutions. They have 
remarkable features which are common to all of them: 

1 . a peak column density of N (H2) 2 x IQ-^^ cm ~^, equivalent 
to Av~ 20 mag (following lBohlin et al.iri978l) . This is low 
compared to most other cases with similar strong C"^0 de- 
pletion but slightly higher than in L1498 whic h has a peak 
column density of N(H2) ~ 1.3 x 10^^ cm"^ (iShirlev et al.l 
'200513 

2. a strong central depletion (^ 30) of C'**0 . 

3. C'**0 depletion starting at a density between 1.5 and 2x10^ 
cm""*. 

4. for both the C'^'O lines (towards the positive offsets) and the 
N2H^ line (for all offsets), the linewidth is so narrow that the 
turbulence contribution to the width is negligible and is set 
to Vtuib(FWHM) ^ 68 ms"'. This is the lowest turbulence 
reported for C'**0 lines to our knowledge. The line width is 
mostly dominated by the thermal width and the enlargement 
due to macroscopic kinematics (rotation and infall). 

5. the radial velocity is strongly constrained in the C'^O enve- 
lope but not so much in the N2H^ internal part. However, a 
constant infall around 0. 1 km s"' is compatible with the data. 
No infall is not. 

6. the rotational velocity is decreasing with increasing radius. 
Solid rotation is therefore excluded and velocity depen- 
dences of the type V oc r^" with a 0.5 gave good solutions. 

7. the temperature pattern requests a slightly higher tempera- 
ture in the N2H^ inner core (~10 K) than in the C"*0 outer 
core (~8 K). If corr ect this would be similar to the B68 case 
jBergin et al.ll2006l) . 

The C'^O temperature is determined by the necessity to fit 
both the (J: 1-0) and (J:2-l) lines with a relatively constrained 
dust profile while the N2H^temperature is the minimum we 
could get trying to minimize also the needed density. Despite the 
observation of a single transition for N2H^, we are constrained 
by the hyperfine structure (Pagani etal. 2007, see ). However 
this temperature difference could possibly prove to be inaccurate 
if the actual H2-N2H^ collisional coefficients were available: 
they are expected to be higher than the present He-N2H^ ones 
and therefore would probably allow to fit the observations with 
both a lower temperature and a lower densi t y. In a ny case, the 
gas appears cooler than what IStepnik et al.l (l2003h have found 
for the dust (~ 12 K). Either their resolution is too coarse and 
mix up a large fraction of the (warmer) envelope emission with 
the core emission, or their single temperature fit is too simple 
by lack of constraints (we have shown that in the case of L183 
cold core the 200 yum emission is not tracing the same dust as 
the s ubmm emissi on and that a temperature gradient is compul- 
sory, l^ganreTal] [2004) or, finally the gas and the dust being 

Note that lShirlev et alj l |2005l) disagree on the total dust amo unt and 
peak density of L1498 with respect to what lTafalla et alj j2004) find. If 
we follow the latter, L1498 is closer to a normal prestellar core than 
to a simple core. T he discrepancy has only partly been explained by 
IShirlev et all ( 120051) . 
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Fig. 9. Fit of the individual 
C"*0 (J: 1-0) spectra along the 
PRONAOS cut. The histogram 
plot represents the data, the 
thick continuous curve, the 
model. A FWHM turbulent 
velocity of 0.11 kms"' has 
been used all over the cut which 
explains why the model lines 
are larger than the observations 
for the positive offsets. Offsets 
along the PRONAOS cut in arc- 
seconds are given in the upper 
right corner of each spectrum. 
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Fig. 10. Same as Fig.|9]for C''^0 
(J:2-l) spectra. 



decoupled at these low densities (n(H2) < 5 x 10"* cm"^ or less), 
the dust could be warmer than the gas though this is not usually 
expected. 

The modeling of the N2D^ (J: 1-0) line indicates a column 
density of 4 to 8 x 10" cm"^, depending on the chosen N2D^ 
abundance profile which we cannot constrain with only one ob- 
servation. The peak N2H^ column density is ~2 x lO'"' cm"^ 
and ther efore the deuterium enric hment is 2 - 4 %. Compared to 
LI 83 (Pagani et al.ll2007l l2009b). we find a similar deuteration 
level at similar densities indicating that the evolutionary speed of 
the deuterium enrichment is probably similar and therefore can 
be used to determine the age of the core. 



Beyond the core characterized by C'^O narrow lines, the en- 
velope, as traced by '^CO, is best described with an almost con- 
stant density (in the range 200-1000 cm"-') and a large turbu- 
lence (Vturb(FWHM) - 0.4 kms '). To follow its macroscopic 
velocity drift we had to introduce a counter-rotating pattern in 
the external layers set constant to -0.3 kms"'. We do not claim 
that the '''CO is indeed rotating but this is the only way to ap- 
proach the velocity profile of this line with our model. Similarly, 
to reproduce the redshifted peak, we had to introduce an expan- 
sion velocity of 90 m s"', that is comparable in amplitude to the 
infall velocity of the core. The fit is not perfect inasmuch as we 
could not centre the peak emission in our model at the same po- 
sition as in the observed position-velocity diagram (Fig. |2ll but 
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Fig. 11. Fit of the individ- 
ual '^CO (J: 1-0) spectra along 
the PRONAOS cut. The his- 
togram plot represents the data, 
the thick continuous curve, the 
model. 
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Fig. 12. Fit of the individ- 
ual NaH^ (J: 1-0) spectra along 
the PRONAOS cut. The his- 
togram plot represents the data, 
the thick continuous curve, the 
model. Only the 3 central hyper- 
fine structure components are 
displayed. A turbulent velocity 
of 0.068 km s ' has been used. 



the fit bears some ressemblance with the reality as shown also in 
Fig.[TT] The derived parameters are therefore reasonably close to 
the actual situation and this envelope is clearly decoupled from 
the core in terms of density profile as much as in terms of turbu- 
lence and large scale motions. That region is also marked by a 
very high abundance of '''CO (6 x 10"^) and the absence of C"*0 
(< 5 X 10"', an abundance of 1 x 10"^ is marginally consistant 
with the observations). The density is not well constrained in 
the external envelope but as the '^CO abundance is already high 
(standard value is 1-2 x 10"^, i.e. 3 to 6 times less), the density 
cannot be lower than what we have assumed here as a lower den- 
sity would increase even more the '^CO abundance. Conversely, 
we could increase the density to lower the abundance of '^CO, 
but it would become higher than the density in the last C"^0 
layers. In the absence of other '^CO transitions it is difficult to 
constrain correctly the density and the abundance in these last, 
turbulent layers. 



The collapsing core has a size of ~3 x 10 AU in radius 
and holds a mass of ~4 Mq. This size is larger than the one 
of L1498 which was already shown to be the largest by 50% 
amo ng the nearby prestellar cores mentioned in the study by 
Shi rley et alj (|2005|) . Our mass estimate is somewhat approxi- 
mate as the density of the N2H^ region is possibly overestimated 
but on the other hand, the core could be slightly more extended 
than a sphere, closer to a cylindar. The estimate derived from 
the dust observations is similar but less constrained due to the 
unknowns on the dust emissivity and temperature in the core 
at scales smaller than the PRONAOS resolution. The '^'CO low 
density filament holds at least 20 Mq accounting only for the 
part along its main axis which covers the core. Its diameter is 
250,000 A.U. or 1.25 pc. 
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Table 2. The Monte-Carlo cloud model used to fit the spectra (Figs.|7] andl9]-fT2li 



Shell Radius 


Density 


Temperature 


Rad. velocity" 


Rot. velocity 


Turbulence'' 


Abundance^ 




(A.U.) 


(") 


(cm"') 


(K) 


(kms-') 


(kms-') 


(km s-') 


N2H+ 


C'»0 


"CO 


1740 


12 


50000 


10 


0.11 


0.200 


<0.11/<0.068 


1.5(-11) 


<5(-9) 


<3(-8) 


3480 


25 


47400 


10 


0.11 


0.141 


<0.11/<0.068 


2.5(-10) 


< 5(-9) 


<3(-8) 


5220 


37 


41300 


10 


0.11 


0.115 


<0.11/<0.068 


4(-10) 


< 5(-9) 


<3(-8) 


6933 


50 


38100 


10 


0.11 


0.100 


<0.11/<0.068 


5(-10) 


<5(-9) 


<3(-8) 


8666 


62 


33900 


10 


0.11 


0.089 


<0.11/<0.068 


5(-10) 


<5(-9) 


<3(-8) 


10400 


75 


24300 


10 


0.11 


0.082 


<0.11/<0.068 


4(-10) 


< 5(-9) 


<3(-8) 


12133 


87 


18500 


10 


0.11 


0.076 


<0.11/<0.068 


3(-10) 


l(-8) 


6(-8) 


13933 


100 


14700 


10 


0.11 


0.070 


<0.11/<0.068 


2(-10) 


1.7(-7) 


l(-6) 


15666 


112 


11900 


9 


0.11 


0.066 


<0.11/<0.068 


I(-IO) 


1.7(-7) 


l(-6) 


17400 


125 


9990 


8 


0.11 


0.063 


<0.11/<0.068 


I(-IO) 


1.7(-7) 


l(-6) 


19133 


137 


8540 


8 


0.11 


0.060 


<0.11/<0.068 




1.7(-7) 


l(-6) 


20866 


150 


7400 


8 


0.11 


0.058 


<0.11/<0.068 




1.7(-7) 


l(-6) 


22600 


162 


6100 


8 


0.11 


0.055 


<0.11/<0.068 




1.7(-7) 


l(-6) 


24333 


175 


4800 


8 


0.11 


0.053 


<0.11/<0.068 




1.7(-7) 


l(-6) 


26066 


187 


3850 


8 


0.11 


0.052 


<0.11/<0.068 




1.7(-7) 


l(-6) 


27866 


200 


3130 


8 


0.11 


0.050 


<0.11 /<0.068 




1.7(-7) 


l(-6) 


29600 


212 


2560 


8 


0.11 


0.048 


<0.11/<0.068 




1.7(-7) 


l(-6) 


31333 


225 


2140 


8 





-0.3 


0.4 




<5(-9) 


6(-6) 


46666 


335 


1000 


8 


-0.09 


-0.3 


0.4 




<5(-9) 


6(-6) 


66666 


479 


500 


10 


-0.09 


-0.3 


0.4 




<5(-9) 


6(-6) 


83333 


598 


400 


12 


-0.09 


-0.3 


0.4 




<5(-9) 


6(-6) 


100000 


718 


400 


13 


-0.09 


-0.3 


0.4 




<5(-9) 


6(-6) 


126666 


909 


200 


13 


-0.09 


-0.3 


0.4 




<l(-9) 


l(-8) 



" positive velocity indicates inward motion (contraction). 
FWHM turbulence is <0.068l£m s ' for the positive offsets in the C '**0 model and for all offsets in the NiH^ model and < 0. 1 1 km s"' otherwise 
(for Figs.|9]and[T0]only one value could be used: 0. 1 1 km s"' has been selected), 
relative to H2. I(-IO) means 1 x 10"'". 
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Fig. 13. H2 density and molecular abundance (relative to H2) as 
a function of radius. 

5. Discussion 

Several aspects of this part of the L1506 filament are sur- 
prising and presently difficult to explain. The expanding and 
counter-drifting external envelope is one of them. One partially 
similar case has however already been reported: Young et al.l 
(Il981h have shown that the B5 core as traced by C'^O had 
a position/velocity profile gradient opposite to the '^CO one 
which they interpreted as a counter-rotating core (but they don't 
discuss the radial c omponent). Contemporary theoret ical work 
dGilhs et al.1 Il979l: iMouschovias & Paleologoul Il979h was in- 



voked to explain this situation: a radial, frozen-in magnetic field 
perpendicular to the rotation axis could brake down the rotation 
of the cloud but because the field remains anchored to the exter- 
nal regions, the magnetic torques do not vanish immediately and 
counter-rotation is induced in the core which could oscillate sev- 
eral times before coming to a halt. Whether this is a similar case 
here is not yet clear because 1) we are not certain that the '^CO 
itself is tracing a rotation, 2) the expansion of the '^CO envelope 
is not predicted by the model and, 3) we know nothing about the 
magnetic field inside this core (external parts have been mapped 
though, Moneti et al.. 1984; Heiles 2000) . 

The strong depletion of C'^O with its abrupt appearance at 
relatively low density is another one. It is similar to the L1498 
case but in sheer contrast with other very young cores, like 
L1495B, L1521B, and L1521E which show no depletion despite 
a density which is typically 5 times higher (and possibly 10 times 
if He-N2H^ coUisional coefficients are indeed over-evaluating 
the actual density in L1506C or if we c ompare to L1498 peak 
density as derived bv' Shirlev et al.ll2"005h . Based mainly on this 
absence of depletion, these cores have been judged to be young 
but it is somewhat strange that having very recently acquired 
their present density state, they show no kinematical sign of 
contraction. It is thus conceivable that they are older than as- 
sumed and that the absence of depletion could be due to another 
mechanism , like a strong desorption, an ill-understood mecha- 
nism dRoberts et alj|2007b . Indeed, it is thought that energetic 
cosmic ray i mpacts can de sorb the most volatile ices like CO 
(lLegeretal.|[r9 85; Hasegaw a & Herbstlll993l:lBringa & JohnsonI 
2004; Roberts et al.i2007l) . The desorption efficiency is inversely 
correlated to the size of the grains as smaller grains are easier 
to warm up. Therefore the question arises whether cores like 
L1495B have smaller grains than L1506C or not and how much 
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different in size should these grains be to allow or prevent an ef - 
ficient desorption of CO. As explained by [Roberts et alj (l2007h . 
all the desorption mechanisms are yet poorly known and it is 
presently difficult to go further. However, if this explanation 
proves to be valid, then one has to understand why grains co- 
agulate in some sources and not in others. Turbulence could be 
one of the main key parameters involved in this differentiation 
between cores. From Chok s hi et a l. ( 1993) and Poppe & Blum 
(Il997h works, iFlower eta show that dust coagulation 

efficiency peaks for turbulent velocities in the range 10-80 
ms"' (see their Fig. A.l), a condition that is hardly met even in 
the most quiescent prestellar cores but fits the present case. On 
the opposite, L1521E shows no depletion while the turbulence 
is suprathermal (even for N2H+, Vtuib(FWHM) = 0.27 kms 
iTafalla&Santiagol 12004 1 a factor of 4 at least larger than in 
L1506C). If the size of the grains in L1521E is not known, large 
grains have b een advocated to e xplain sub-millimeter emission 
in L1506C: Stepnik et all (l2003b have built a model where the 
very small grains, emitting at 60 jum, disappear for offsets lower 
than ro = 4' from the centre while the sub-mm emissivity in- 
creases by a factor of 3.4 due to the formation of fluffy ag- 
gregates. Their observations were taken with an angular reso- 
lution of 2-3.5', so the value of ro could be lower and fit the 
extent (+ 200") of the low turbulence region in L1506C. Grain 
growth, depletion and desorption are certainly not completely 
understood at the moment. 

The absence of C'^O on 5/6th of the external envelope width 
as traced by the POM-1 '^'CO observations is another puzzle. 
The '^CO abundance itself seems to be higher than usual so 
that the '^CO/C"^0 ratio reaches 600-1200 over a large region. 
Though '^CO fractionation does exist, it is weak and does not 
seem to be able to explain this difference. More observations of 
the envelope (including '^CO (J: 1-0) and (J:2-l) and very low 
noise C'**0 transitions) are needed to address this issue. 

Clear trace of collapse among prest ellar cores has been found 
in a few cases. iLee et al.] (I1999L 12004 ) observed a large num- 
ber of asymmetrically blue profiles towards a sample of starless 
cores but they traced the CS kinematics with respect to N2H^ 
or DCO^ and, as CS is depleted inside the N2H^ region, they 
trace the inward motion of the envelope only. A few prestellar 



one being LI 544 (e.g. 


WilUams et al.lll999l: iTafalla et alJl2002l 


Ivan der Tak et al.ll2005 


) which was reported to clearly undergo 



voke the possibility of infall from the observations of two other 
cores, L1498 and L1517B despite t he absence of a clear velocity 
gradient and Willi ams et al.l (l2006h found that L694-2 has a pro- 
file reminiscent of LI 5 44, indica ting also probable infall. This 
has been confirmed by iLee et al., (200?) who also report infall 
of LI 197. Here the inward motion of N2H^ is not strongly con- 
strained in amplitude but unavoidable to reproduce the velocity- 
profile plot: no infall at all gives clearly a wrong result but the 
details of the radial velocity profile are uncertain. The inward 
motion is also clearly shown by the C"^0 andC'^O lines. It is the 
first time to our knowledge that such lines do trace the collapse 
of a core . We think that this is possible here because several 
favorable factors are met : 

1. no CO isotopologues in the central part, separating clearly 
the front and rear moving layers as indicated by C'^O . 

2. no C'^O in the extended envelope which would have hidden 
or blurred the faint variation of intensity between the blue 
and red emissions which let us differentiate between expan- 
sion and contraction 



3. low C'**0 velocity dispersion, increasing the line opacity 
near unity (ry i-o = 0.56, tj-2-i = 0.91). Combined with 
the outward decreasing excitation temperature, it allows for 
the weak but measurable differential absorption between the 
blue and red components. 

The most surprising aspect is that this core with a density 
below ~1 X 10^ cm is supposed to be dynamically stable 
given its total mass (less than 10 Mm. iKeto & Casellil l2008h 
which is clearly not the case here (the core mass being es- 
timated to be ~4 Mq). However, as noted by iKeto & Casellil 
(2008), most clouds without non-thermal energy support would 
be unstable and in this case, turbulence is abnormally low, non- 
measurable in several parts of the core. This would support the 
idea of a possible collapse though we have no clue on the mag- 
netic field intensity and direction. If a magnetic field is indeed 
present to induce the core counter-rotation, it is not incompati- 
ble with infall as it is included in the magnetic braking model of 
iMouschovias & Paleologoiil (11979 ). Oscill ations or so und waves 
have been invoked in the case of e.g. B68 (Maret et al.>200 7) and 
one could be tempted to invoke a similar argument here but the 
rotation of C'^O and "counter-rotation" of ^^CO seem difficult 
to explain by sound waves. 

6. Conclusions 

The large size of L1506C and its low density clearly indicate 
that it is not yet a prestellar core but its inward motion and its 
kinematical decoupling with the outer envelope are signs that it 
is in the process of turning into one. It is the first time that a low 
density core is clearly observed to collapse to form a prestel- 
lar core, both the kinematical and density status of L1498 being 
somewhat unclear. Two different evolutionary paths from nor- 
mal clouds to prestellar cores seem to exist : one which starts 
with density enhancement first much before any depletion occurs 
and one which starts with depletion first. Two questions arise : 
what is the reason for this differentiation (turbulence via its ac- 
tion on grain growth, itself limiting the desorption capability of 
cosmic rays ?) ? and are both paths really leading to standard 
prestellar cores ? Coming observations with the Herschel Space 
Observatory will help to address these questions and the nature 
of ffie dust inside ffiis object. 

Acknowledgements. We want to thank an anonymous referee for her/his fruit- 
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fruitful discussions. 



References 

Abergel, A., Boulanger, F., Mizuno, A., & Fukui, Y. 1994, ApJ, 423, L59 
Bergin, E. A. & Langer, W. D. 1997, ApJ, 486, 316 
Bergin, E. A., Maret, S., van der Tak, F F S., et al. 2006, ApJ, 645, 369 
Bernard, J. P., Boulanger, F, Desert, F X., & Puget, J. L. 1992, A&A, 263, 258 
Bernes, C. 1979, A&A, 73, 67 

Bohlin, R. C, Savage, B. D., & Drake, J. F 1978, ApJ, 224, 132 
Bringa, E. M. & Johnson, R. E. 2004, ApJ, 603, 159 
Caselh, P 2002, Planet. Space Sci., 50, 1133 

Chokshi, A., Tielens, A. G. G. M., & Hollenbach, D. 1993, ApJ, 407, 806 
Daniel, F, Dubemet, M.-L., Meuwly, M., Cemicharo, J., & Pagani, L. 2005, 

MNRAS, 363, 1083 
Dupac, X. & Giard, M. 2002, MNRAS, 330, 497 
Ehas, J. H. 1978, ApJ, 224, 857 

Flower, D. R., Pineau Des Forets, G., & Walmsley, C. M. 2005, A&A, 436, 933 
Gillis, J., Mestel, L., & Paris, R. B. 1979, MNRAS, 187, 311 
Hasegawa, T. I. & Herbst, E. 1993, MNRAS, 261, 83 
Heiles, C. 2000, AJ, 119, 923 

Hirota, T., Ito, T., & Yamamoto, S. 2002, ApJ, 565, 359 
Hirota, T., Maezawa, H., & Yamamoto, S. 2004, ApJ, 617, 399 



12 



L. Pagani et al.: L1506 : a prestellar core in the making 



Kenyon, S. J., Dobrzycka, D., & Hartmann, L. 1994, AJ, 108, 1872 
Keto, E. & Caselli, P. 2008, ApJ, 683, 238 
Klapper, G., Surin, L., Lewen, R, et al. 2003, ApJ, 582, 262 
Lada, C. J., Huai'd, T. L., Crews, L. J., & Alves, J. F. 2004, ApJ, 610, 303 
Lee, C. W. & Myers, P. C. 1999, ApJS, 123, 233 
Lee, C. W., Myers, P C, & Plume, R. 2004, ApJS, 153, 523 
Lee, C. W., Myers, P C, & Tafalla, M. 1999, ApJ, 526, 788 
Lee, S. H., Pai'k, Y.-S., Sohn, J., Lee, C. W., & Lee, H. M. 2007, ApJ, 660, 1326 
Leger, A., Jura, M., & Omont, A. 1985, A&A, 144, 147 
Maret, S., Bergin, E. A., & Lada, C. J. 2007, ApJ, 670, L25 
Moneti, A., Pipher, J. L., Heifer, H. L., McMillan, R. S., & Perry, M. L. 1984, 
ApJ, 282, 508 

Mouschovias, T. C. & Paleologou, E. V. 1979, ApJ, 230, 204 
Myers, P C, Mardones, D., Tafalla, M., Williams, J. P, & Wilner, D. J. 1996, 
ApJ, 465, L133+ 

Nercessian, E., Castets, A., Benayoun. J. J., & Cemicharo. J. 1988. A&A, 189, 
207 

Onishi, T., Mizuno, A., Kawamura, A., Ogawa, H., & Fukui, Y. 1996, ApJ, 465, 
815 

Ossenkopf, V. & Henning, T. 1994, A&A, 291, 943 
Pagani, L., Bacmann, A., Cabrit, S., & Vastel, C. 2007, A&A, 467, 179 
Pagani, L., Bacmann, A., Motte, R, et al. 2004, A&A, 417, 605 
Pagani, L. & Breart de Boisanger, C. 1996, A&A, 312, 989 
Pagani, L., Daniel, P., & Dubernet, M.-L. 2009a, A&A, 494, 719 
Pagani, L., Pardo, J.-R., Apponi, A. J., Bacmann, A., & Cabrit, S. 2005, A&A. 
429, 181 

Pagani, L., Vastel, C, Hugo, E., et al. 2009b, A&A, 494, 623 
Poppe, T. & Blum, J. 1997, Advances in Space Research, 20, 1595 
Roberts, J. R, Rawlings, J. M. C, Viti, S., & Williams, D. A. 2007, MNRAS, 
382, 733 

Shirley, Y. L., Nordhaus, M. K., Grcevich, J. M., et al. 2005, ApJ, 632, 982 
Stepnik, B., Abergel, A., Bernard, J.-P, et al. 2003, A&A, 398, 551 
Tafalla, M., Myers, P C, Caselh, P, & Walmsley, C. M. 2004, A&A, 416, 191 
Tafalla, M., Myers, P C, Caselli. P, Walmsley, C. M., & Comito, C. 2002, ApJ, 
569, 815 

Tafalla, M. & Santiago, J. 2004, A&A, 414, L53 
van der Tak, R. R S.. Caselli, P, & Ceccarelli, C. 2005, A&A, 439, 195 
Willacy, K., Langer, W. D., & Velusamy, T. 1998. ApJ, 507, L171 
Williams, J. R, Lee, C. W., & Myers, R C. 2006, ApJ, 636, 952 
Wilhams. J. P, Myers, P C. Wilner, D. J., & di Francesco, J. 1999, ApJ, 513, 
L61 

Young, J. S., Goldsmith, P R, Langer, W. D., & Wilson, R. W. 1981, ApJ, 251, 
LSI 

Zucconi, A., Walmsley, C. M., & Galli, D. 2001. A&A, 376, 650 



List of Objects 

'L1506' onpageH 
'L1506C' onpagell 
'B68' on page[I[ 
'L183' on pagey. 
'L1498' on page 
'L1544' on page 



'L1495B' on page 
'L1521B' on page 
'L1521E' on page, 
'L1498' onpage[ir 
'G2' on page[T] 
'L1506' on page|2 
'L1506C' on pagell 
'L183' on pagej?! 
'B5' on page[l 
'L1544' on page 
'L1498' on page 



11 
11 

'L1517B' on pagejll 
'L694-2' on pageJuT 
'LI 197' on page[lir 



